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An issue currently under debate in the literature is how far from the black hole is the Fermi-observed GeV 
emission of powerful blazars emitted. Here we present a clear diagnostic tool for testing whether the GeV 
emission site is located within the sub-pc broad emission line (BLR) region or further out in the few pc scale 
molecular torus (MT) environment. Within the BLR the scattering takes place at the onset of the Klein-Nishina 
regime, causing the electron cooling time to become almost energy independent and as a result, the variation of 
high-energy emission is expected to be achromatic. Contrarily, if the emission site is located outside the BLR, 
the expected GeV variability is energy-dependent and with amplitude increasing with energy. We demonstrate 
this using time-dependent numerical simulations of blazar variability. 


1. Introduction 

Blazars are by far the most common objects de- 
tected in the gamma-ray sky [Abdo et al. 2011]. Fermi 
has detected blazar variability as short as a few hours 
[e.g. Abdo et al. 2010]. During these flares, the GeV 
luminosity has been known to increase by a factor 
of up T.O several compared to its pre-flare luminosity. 
Because blazars cannot be resolved at these energies 
(or at any other energy, with the possible exception 
of VLB A observations), it is impossible to determine 
the location of these flares by direct detection. 

To address this issue, we propose a diagnostic test 
that utilizes Fermi variability data of short flares to 
determine the location of the GeV emission in blazars. 


2. Sources of Seed Photons 

Relativistic effects determine which photon field is 
dominant at varying distances from the central black 
hole, and as a result the location of the GeV flar- 
ing site determines the dominant source of seed pho- 
tons. The co-moving (jet frame) energy density of a 
radiation field U* scales as differing factors of T de- 
pending on the direction from which the photons en- 
ter the emitting region [Dermer & Schlickeiser 1994, 
Georganopoulos et al. 2001]. If the photon field is 
isotropic U' ~ Y 2 U. For photons entering the emit- 
ting region from behind the relativistically moving 
blob U f « UT~ 2 . If we assume a nominal FSRQ ac- 
cretion disk luminosity of L^isk ~ 10 45 erg s -1 and 


that a fraction £ = 0.1 of this radiation is reprocessed 
by both the BLR and the MT, a typical luminosity 
of the external radiation field is L ext ~ 10 44 erg s -1 
[Ghisellini & Tavecchio 2009], 

If the emission site is located within the BLR (at 
R ~ 10 17 cm), the photon field can be considered 
isotropic in the galaxy frame and its co-moving energy 
density is U' BLR ft* 2.6 T^ 0 L blrm^blr , 17 er S cm -3 . 
Similarly, the AIT photon field is isotropic inside the 
BLR and its co-moving seed photon energy density 
is U’ MT ft* 2.6 x 10 -2 T\ 0 L M tm r 'mt,i% er S cm -3 . 
Clearly, inside the Rblr the co-moving BLR photon 
field energy density dominates over that of the AIT by 
a factor of ~ 100. 

If the emission site is located at R ~ 10 18 cm 
(within the MT) then the BLR UV photons enter 
the emitting region practically from behind, so that 
^ BLR w 2.6 X 10 4 L B L R, 4:4 RfrfT, 18^ 10 er S cm 3 - 
The IR photons from the MT retain the same co- 
moving energy density previously given by U' MT & 
2.6 x 10 -2 er S crn_3 * In this case, 

therefore, it is the MT that dominates the co-moving 
photon energy density. 

These external photon field co-moving luminosi- 
ties need to be compared to the synchrotron pho- 
ton field energy density. If Ru 0 b is the size of the 
emitting blob, the co-moving synchrotron photon en- 
ergy density is U$ r 4 • The most P lausi " 

ble assumption for the size of the emitting region, 
however, is to set an upper limit to it by its vari- 
ability timescale: Ruob = ct var 5. We then obtain a 
lower limit for the synchrotron energy density U f s & 
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6.3 x 10 -2 L 46 £“a r6 / l IT 0 6 erg cm" 3 , where we used a 6 
hour variability scale, as seen in variability timescales 
observed by Fermi [Abdo et al. 2010], 

Although the synchrotron photon energy density is 
substantially lower than the BLR photon energy den- 
sity when the blazar emission site is within the BLR, 
it is comparable or even higher than the MT pho- 
ton energy density when the blazar emission site is 
within the MT but outside the BLR. We briefly note 
here if we want to have a blob of a given size Rbiob 
at a distance larger than R/T , [e.g. Marscher et al. 
2008] the blob cannot occupy the entire cross section 
of the jet. Our diagnostic hinges on the fact that if 
the blazar emission site is located outside the BLR, 
electron cooling occurs in the Thomson regime. This 
is always true for the case of the MT photons being 
dominant and it can be shown that if SSC photons 
dominate, cooling still takes place within the Thom- 
son regime for powerful FSRQs. 


3. Cooling in the BLR vs Cooling in the 
MT 

The critical difference between the BLR and the 
MT is the energy of the seed photons : photons orig- 
inating from the BLR axe UV photons (co « 10~ 5 ) 
while photons originating from the MT are IR pho- 
tons (to ~ 10~ 7 ). This difference by a factor of ~ 10 2 
in typical photon energy is critical in that it affects the 
energy regime in which electron cooling takes place , 
and thus the energy dependence of the electron cool- 
ing time. 

In powerful FSRQs the IC emission in high states 
can dominate over the synchrotron emission by a fac- 
tor of up to ~ 100 [e.g. Abdo et al. 2010]. In cases of 
high Compton dominance, the primary electron cool- 
ing mechanism is IC scattering. For electrons cool- 
ing in the Thomson regime (7^0 < 1, where both the 
electron Lorentz factor 7 and the seed photon energy 
e 0 are measured in the same frame) the cooling rate 
7 oc 7“. For electrons with 760 ^ 1 cooling takes place 
in the Klein-Nishina (KN) regime with 7 oc In 7 [Blu- 
menthal &; Gould 1970]. Given that the BLR photons 
are ~ 10 2 more energetic than the MT photons, we 
expect that if cooling takes place in the BLR, effects 
of the transition from the Thomson to the KN regime 
will be manifested in electrons of energy ~ 10 ” 2 times 
lower than if the cooling takes place within the MT. 

The effects of the transition between Thomson and 
KN regimes on the electron energy distribution (EED) 
and the resultant spectrum of the synchrotron and IC 
emission have been studied before [e.g. Blumenthal 
1971, Dermer & Schlickeiser 1993, Georganopoulos et 
al. 2006, Kusunose & Takahara 2005, Manolakou et al. 
2007, Moderski et al. 2005, Sikora et al. 2009, Sokolov 
et al. 2004, Zdziarski 1989]. In short, because 7 oc 




Figure 1 : Cooling time in the blob frame as a function of 
co-moving electron energy 7. Top panel: Blazar emission 
site located in the BLR. Bottom panel: Blazar emission 
site located in the MT. The dotted lines represent the 
various cooling mechanisms (blue = Thomson cooling, 
red = KN cooling, green = synchrotron cooling), the 
solid black line is the total cooling time. Plots were 
calculated for the following values: seed photon energies: 
to, BLR = 3 x 10” 6 and co.mt = 6 x 10 -7 , energy 
densities Ublr = 2.65 x 10” 2 erg cm -3 , Umt = 3x 10 -4 
erg cm” 3 , a Ub that corresponded to a Compton 
dominance of ~ 100, and Tbuik = 20. 

7 2 in the Thomson regime and 7 oc In 7 in the KN 
regime, the cooling time Tcooi = 7/7 scales as 7 -1 in 
the Thomson regime and as 7/ In 7 in the KN regime. 

Because the cooling time is approximately energy 
independent around 760 ~ 1, this energy-independent 
cooling time will be manifested at energies lower by 
a factor of ~ 100 for cooling taking place in the BLR 
compared to cooling taking place in the MT, since the 
BLR seed photons have an energy higher than that of 
the MT by a factor of ~ 100. This can be seen in Fig. 
1 where we plot the electron cooling time for a source 
with a ratio of external photon field energy density Uq 
in the co-moving frame to co-moving magnetic field 
energy density Ub, Uq/Ub = 100. This corresponds 
to a factor of ^ 100 Compton dominance (ratio of 
inverse Compton to Synchrotron luminosity), similar 
to what is observed in the most Compton dominated 
sources. 

The transition from Thomson cooling to KN cool- 
ing, and from KN cooling to synchrotron cooling can 
also be seen in Fig. 1. If the electrons are cooling 
on photons from the BLR, cooling takes place in the 
Klein-Nishina regime and the cooling time scale is ap- 
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proximately energy independent around 7eo rw 1 (Fig. 
1, top panel). If the electron population cools on pho- 
tons from the molecular torus, cooling takes place in 
the Thomson regime (Fig. 1, bottom panel). The 
cooling time is heavily energy dependent, and any 
variations should consequently exhibit heavy energy- 
dependence. 


4. The Diagnostic Test 

The energy dependence of the cooling time results 
in an energy dependence (or lack thereof) of varia- 
tions: if the blazar emission site is located within the 
BLR, variations should be achromatic. The energy 
dependence of the variations can be used as a diag- 
nostic test to determine if the GeV emission site is 
located within the BLR. By comparing Fermi light 
curves of flares at different energies, we propose that 
the energy dependence of the light curve can be used 
as a diagnostic test to rule out whether the GeV flare 
originates in the BLR or MT. 


4.1. Numerical Simulation Results 

To demonstrate the effect of the energy indepen- 
dence or dependence of the electron cooling time on 
the variability of a flare, we utilized a one-zone nu- 
merical model to simulate a flare. We initialized the 
code with values appropriate for a high power blazar 
with s Compton dominance of ~ 100. For this partic- 
ular simulation we assumed a source size R = 10 16 cm, 
bulk Lorentz factor T = 10, co-moving injected elec- 
tron luminosity L = 2 x 10 44 erg s -1 , maximum 
electron Lorentz factor 7 mar , = 10 5 , and electron in- 
dex p = 2.5. For the case of a flaring region located 
within the BLR we assumed an initial photon energy 
€0 = 3 x 10 -5 and an energy density (in the galaxy 
frame) Ublr = 2.6 x 10” 2 erg cm -3 . For the case of 
a flaring region located outside the BLR we assumed 
an initial photon energy eo = 6 x 10~ 7 and an energy 
density Umt = 3 x 10 -4 erg cm -3 . For each case 
the magnetic field B was fixed to assume a Compton 
dominance Uec/Ub = 100. 

The system behaves as expected, showing a notice- 
able difference in the decay rate as well as the am- 
plitude of the flare depending on if the seed photons 
originated from within or from outside the BLR (see 
Figs. 2 and 3). This difference in decay rate and 
amplitude can be used as a diagnostic test to differen- 
tiate between flares that take place inside or outside 
the BLR. 

4.2. Light-Crossing Time Effects 

Because the GeV emitting region is not a point 
source, any change in the luminosity of the blob will 



Figure 2: Light curve at various energies (BLR seed 
photons): e 0 = 100 MeV, eo = 1 GeV , eo = 10 GeV 



Figure 3: Light curve at various energies (MT seed 
photons): eo = 100 MeV, eo = 1 GeV, eo — 10 GeV. 


not be seen instantaneously. Instead, the observed de- 
cay time of the light curve is the result of the actual 
decay time and the light-crossing time inherent in the 
blob. To test whether the light-crossing time would 
erase any difference in the light curves (as predicted 
by our diagnostic) for the case of the emitting region 
being located outside the BLR, we modeled the flux 
of a source with light-crossing time that is de- 
creasing in flux. We assumed the flux of the source is 
decaying exponentially, F(t) = F 0 exp -t / tc , where t c 
is the cooling time at a specific energy. 

Our diagnostic predicts differences in the decay time 
of the light curves for the case where the emitting 
region is located within the MT; for the purposes of 
this demonstration we assume cooling occurs in the 
Thomson regime and as a result t c oc e -1 / 2 . We plot 
the resultant light curves for cooling times differing by 
a factor of l/\/l0 (i.e. energies differing by a factor 
of 10). As evident in Fig. 4, even with light-travel 
time effects convolved with exponential decay times, 
the predicted differences in the decay times are still 
preserved. 
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Figure 4: Light curves with light-crossing time effects 
factored in. Initial flux and light-crossing time are 
normalized to 1. Light curves are plotted for energies 
differing by a factor of 10. 

4.3. Feasibility Study 

A flare that occurs outside the BLR should exhibit 
energy dependent cooling times. Practical applica- 
tion of the diagnostic hinges on the requirement that 
the decay times in different energy bands should have 
detectable differences. To demonstrate this, we took 
as a test case a flare of 3C 454.3 [Abdo et &L 2011] 
and assume an exponential decay of the light-curve 
F(e, t) = Nf, -f iVoe -< / <c , where N & and Nq are nor- 
malizing factors and t c is the energy-dependent cool- 
ing time. We assumed that the flaring site is located 
outside the BLR (i.e. cooling occurs in the Thom- 
son regime) because this is case where different decay 
times need to be resolved. We created simulated data 
points at 6-hour time intervals, applied the energy- 
dependent maximum error from the Fermi observa- 
tions of the flare, and then fit a curve to the simu- 
lated data to obtain a cooling time. From this test, we 
found that even within maximum observational error, 
the decay times at different energies are still distin- 
guishable. 


5. Conclusions 

We have presented a diagnostic test that utilizes 
blazar variability to determine the location of the GeV 
emitting site in blazars. The energy difference in seed 
photons originating from the BLR versus seed photons 
originating from the MT causes electrons within the 
emitting site to cool in different energy regimes. 

For the case where the GeV emitting site is located 
within the BLR, cooling takes place at the onset of the 
KN regime, and the resultant electron cooling time 
is energy-independent. We have demonstrated that 
the associated light curves exhibits decay times that 
are approximately energy independent. Conversely, 
for the case where the GeV emitting site is located 
outside the BLR, cooling takes place in the Thom- 
son regime and the electron cooling times are heavily 


energy dependent. In this case, the associated light 
curves exhibit energy dependence of their decay times. 

The energy dependence of the decay time of the 
light curves is visible within the Fermi energies; these 
differences can be used as a diagnostic test to deter- 
mine whether the GeV emitting region is located in- 
side or outside the BLR. These effects are observable 
within the maximum measured error of Fermi obser- 
vations and are not erased due to considering light- 
travel time effects. If light curves from a sufficiently 
bright and rapid flare [such as that in 3C 454.4; Abdo 
et al. 2011] are compared at different energies, if the 
GeV emitting site is located within the BLR, the de- 
cay times will exhibit no energy dependence, whereas 
if the emitting site is located within the MT, the decay 
times will exhibit energy dependence. 
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